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Solar energy has become one of the most studied topic in the field of 
renewable energy. In this paper, an artificial intelligent approach is proposed 
for the optimization of a photovoltaic solar energy harvesting system. An 
Electromagnetism-Like Mechanism Algorithm (EM) has been developed to 
search for the hourly optimum tilt angles for photovoltaic panels. In order to 
investigate the effect of the search step size on the efficiency and overall 
accuracy of the algorithm, the EM has also been modified into several 
variants with different search step size settings. Experimental findings show 
that EM with bigger search lengths has the advantage of reaching a near 
optimum tilt angle in earlier iterations but less accurate. EM with smaller 


: step lengths, on the other hand, can hit a relatively more optimum tilt angle 
Photovoltaic panel in the process. During the peak of the power generation at noon, EM with 
Solar tracking system smaller search stes found an optimum tilt angle which yielded additional 
Tilt angle 3.17W of power compared to a fixed panel. We thus conclude that the 
proposed EM performs well in optimizing the tilt angle of a photovoltaic 
solar energy harvesting system. 
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1. INTRODUCTION 

Renewable energy has become a term commonly known worldwide nowadays. With the fossil fuel 
diminishing, researchers around the globe are turning their attention into more sustainable energy sources. 
Among others is the solar energy [1]. Solar energy is popularly researched and implemented around the 
world due to its sustainability and minimal environmental damage [2]. Two systems are widely used in 
harvesting solar power today, namely the Solar Thermal Collector (STC) and the Photovoltaic (PV) cells. 
STCs convert sunlight into thermal energy while the PV cells turn solar energy directly into electricity [3]. 

In the literature, there are many methods of improvement to increase the energy generated by a solar 
harvesting system. One of the important aspects is to maximize the solar radiation that falls onto the 
harvesting system. The solar radiation data are not readily available for many locations in the world. Thus, it 
has to be estimated. A number of techniques can be found in the literature to estimation the solar radiation on 
tilted surface [4]-[6]. The total solar radiation that falls onto a PV panel is greatly affected by its orientation 
and tilt angle [7], [ 8]. Thus, one of the best ways to make sure maximum solar radiation onto the PV panel is 
to implement a solar tracker. 

The contribution of this paper is threefold and can be summarized along the line as follows. 
Firstly, a relatively new approach, the Electromagnetism-like Mechanism algorithm (EM), has been 
implemented to track the maximum solar ray in order to optimize the power generated by a PV solar energy 
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harvesting system. To date, the implementations of the EM algorithm in solar tracking systems remain rather 
limited. Secondly, an analytical study on the effect of the search step length in EM has been carried out to 
investigate the gravity of it to the outcome of the optimization process. Finally, the performance of the 
modified EMs are evaluated and benchmarked. The outline of this paper can be divided into 5 main sections. 
The second section of the paper introduces the trajectory of the solar tracking system used in this research 
and the formulation of the objective function. In Section 3, the EM algorithm is discussed in details. 
Section 4 shows the experimental results, analysis and detailed discussion. The final section offers the 
conclusion. 


2. TRAJECTORY OF THE SOLAR TRACKING SYSTEM 

It is important to make sure a solar panel operates at an optimum tilt angle in order to expose the 
panel as close to its maximum solar intensity as possible. Incorrect positioning can lead to loss of potential 
solar power. In this research, a solar energy harvesting system with solar tracking feature is developed. 
A PV panel with 1.651m length and 0.9906m width is employed. A permanent magnet DC motor is installed 
as the main driving actuator to adjust the tilt angle (@;) of the solar panel within the range of 0° to 90°. 
The system is developed to harvest solar energy in the city of Petaling Jaya, Malaysia (3.0833° N, 101.6500° 
E). An optimization procedure is proposed in this study to optimize the output power P with respect to the tilt 
angle of the PV panel. 

Generally, the total solar radiation on a tilted surface (H7) is the sum of the beam solar radiation 
(Hz), diffuse radiation (Hp), and ground reflected radiation (Hp) onto the tilted surface. Thus, for a surface 
tilted at an angle from the horizontal, the incident total radiation is given by the relation in Equation (1): 


Hr=Hg+Hp+ Hp (1) 


The beam radiation onto the tilted surface is as expressed in Equation (2), where H, and Hy, are the 
monthly mean daily global and diffuse radiations on a horizontal surface. 


Hz = (H.-H) Rp (2) 


Since Malaysia is located on the northern hemisphere, the R, for the surface sloped toward the 
equator is given by Equation (3), in which ¢, 6, and a denote the latitude, the declination angle, and the sun 
hour angle for the tilted surface respectively [9]. 


Ri, = $090=85) cos 6 sinact (1/180) asin(O— 44) sind 2 


cos cos 6 sin a+ (1/180) asinOsind 
The sky-diffuse radiation can calculated using Equation (4): 
Hp = RaHa (4) 
where according to [10], the calculation of Ris as shown in Equation (5). 
Ra= (1 - 8s) / 180 (5) 


Thus, the total solar radiation on a tilted surface is as shown in Equation (6), where p refers to the 
ground albedo [9]. 


1-cos 05 


Hy = (H, — Ha)Ry + Hp = HaRa (6) 


From this, the total solar radiation falling on tilted surface is computed by varying tilt angle (0s) 
from 0° to 90°. The optimum tilt angle is taken at which solar radiation on the tilted surface H; becomes 
maximum. In this research, the adjusting steps of the tilt angle is described below in Equation (7), where k 
refers to the iteration index and y(k) denotes the step size of position angle increment or decrement at k. 


OP (k) 


8s (k +1) = 65 (1) + nd) EO 





(7) 
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It is known from [11] that the output power, P of the PV panel is directly proportional to the 
intensity of the solar radiation (H7), given by Equation (8), where A refers to the surface area of the PV panel. 


P= H,A (8) 

Since the beam solar radiation is the main component in the incident total radiation, the intensity of 
solar radiation can be expressed as Equation (9) [12] where H,, is the maximum solar radiation intensity 
possible, and 6; is the angle between the normal to the panel surface and the sun’s rays. 

Hr = H,,cos(@;) (9) 

It can be observed from Figure | that at any point in time, Equation (10) exists. 

6,+ a+ 8, = 90° (10) 


When the PV panel is tracking the sun, 6;= 0, and correspondingly, a + 4; = 90°. From Equations 
(8), (9), and (10), we know that 


P = H,A * cos(6;) = H,A * sin(6,) (11) 

where 6, = a@ + @s. According to Al-Mohamad (2004), the corresponding solar radiation intensity 
H,, in Equation (9), for a particular solar hour angle a, can be calculated from Equation (12) where a refers to 
the elevation angle of the sun which is varying. 

H, = B * exp[— C/sin(a@)] (12) 

B and C are site- and climate-related constants and are calculated by Equation (13) and Equation 
(14) using data available from the meteorological department and the global constant’s tables published by 
ASHRAE [13], where N is the day of the year and F, is the cloudiness coefficient. 

B=0.132 + 0.023 cos(F.N) (13) 


C = 0.047 + 0.03 cos(F.N) (14) 


Vertical 
' 


Normal 





Figure |. Illustration of angles definition 


3. ELECTROMAGNETISM-LIKE MECHANISM ALGORITHM 

There is a rich literature on the implementation of artificial intelligent algorithms and global 
optimization techniques in solar energy harvesting systems. Optimization algorithms are generally developed 
with the aim to search for the global optima parameter that yields the best optima values in a particular 
optimization problem. The Electromagnetic-like Mechanism (EM) algorithm is a population based global 
optimization search mechanism proposed by Birbil and Fang [14]. Guided by the electromagnetism theory, 
the EM imitates the attraction and repulsion mechanisms of electromagnetic charges in search for the global 
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optimal solution in bounded variables. In the algorithm, each of the charged particles in play is tabulated in 
the possible solution range. The charge magnitude of each particle relates to the objective function value. 
Particles with better objective yields will apply attracting forces while particles with worse objective values 
will apply repulsion forces onto other particles [15]. Bigger difference in objective values generates higher 
magnitude of attraction or repulsion force between the particles. The particles are then moved based on 
superposition theorem. In this research, the EM is applied to maximize the output power P generated with 
respect to the tilt angle, 0,. Figure 2 shows an example of the total force, F, applied on Q, by the repulsive 
force from Q, and attractive force from Q,. 





Figure 2. Total force exerted on Qa by Qb and Qc 


3.1. Conventional EM 
The general flow of a conventional EM is as shown in Table 1. There are five major steps in the 
EM, namely initialization, local search, charge calculation, force calculation, and particles displacements. 


Table 1. General EM Flow 
Algorithm 1: EM (m, MAXITER, LSITER, 6) 
m= number of initial particles 
MAXITER: maximum number of iterations 
LSITER: maximum number of local search 
iterations 








: local search parameter, 6 €(0,1) 
: Initialize () 

: iteration | 

: while iteration < MAXITER do 
Local (LSITER, 6) 

F <— CalcF () 

Move (F) 

iteration < iteration + 1 
: end while 








3.2. Initialization 

In the initialization stage of EM, the feasible ranges of all the tuning parameters (upper bound, 
u, and lower bound, /;,) are defined. Then, a number of initial particles (m) are randomly sampled from the 
feasible solution dimensions, each are taken as an N dimensional hyper-solid. Each value of a dimension in a 
particle is assumed to be uniformly distributed inside the upper and lower bound [16]. In this research, the tilt 
angle is set to be the tuning parameter which varies in the range of 0° to 90°. Thus, /, is set to be 0 while u, is 
set to be 90. The particles are to be evaluated with Equation (11). The particle with the largest objective value 
is marked as the best particle as this is a maximization problem. 


3.3. Local Search 

In this stage, the particles gather local information in the neighborhood and make comparisons. 
The original local search procedure of a conventional EM uses a line search with random search step sizes 
within the feasible range of a solution. This simple line search performed by tunning a particle along its 
dimensions one after another, restricted by a maximum and minimum feasible random step length of A € 
(0,1) [17]. A new random step size is generated for every new iteration. In a conventional EM, this loop is 
immediately exited upon hitting any better objective value. 
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3.4. Charge Calculation 

In this stage, the charge of each and every particle is calculated. This will then follow by the force 
calculation, in which a total force vector exerted onto a particle is obtained using the Coulomb’s Law (Lee et 
al. 2012). The charge of each particle is heavily dependent on its current objective value compared to the best 
solution found. A particle will determine if it exerts attration or repulsion force onto another particle when 
the values of the two charges are compared. The calculation of the charge (q') is shown in Equation (15) 


i ~ £(x')—F (best) ’ 
q = exp nso) et 





where n refers to the number of dimensions in the particle and m represents the size of the 


best 


population. f(x") denotes the best objective value obtained so far. 


3.5. Force Calculation 

The forces generated by one particle onto another can be computed based on the calculated charges 
of each particle. A particle with a relatively better objective value will exert attraction force onto another 
particle while the particle with worse objective value will repulse the other particle. Based on the 
electromagnetic theory, the force of one particle onto another is inversely proportional to the square of the 
distance between the two particles and directly proportional to the product of their charges [18]. The total 
force vector for a particle can be determined by condering the collective forces generated using 
Equation (16). 


a0 fe igi j i 

. (xJ-xé)_ 4 _ if f(x/)<f(x') 
Fi= ji a ae E ; Wi oe 

Gece if f(x)ef(x') 


where f(x’) < f(x') denotes attraction and f(x’) > f(«') refers to repulsion. 


3.6. Particle Movement 

In this stage, all the particles but the best are mobilized to a new location in the feasible solution 
space. This step is crucial to ensure a global exploration of all possible solutions. The the movement of a 
particle is calculated based on Equation (17), where A denotes the global particle movement step length. 
In this research, it is set to be a random value between 0 and 1, assumed to be uniformly distributed between 
the upper boundary (u;,= 90°) and the lower boundary (/,= 0°). 


Xe ob APE (tig = xt) FE >0 
Xp — Xe + AFL (xk - Le) 5 Fe <0 (17) 


Holding the absolute power of attraction towards all other particles, the best particle does not 
displace [19]. After a pre-determined number of iterations, the best tilt angle found by the EM is then fed into 
the actuator control to tilt the solar panel accordingly. 


3.7. The modified EM 

The setting of the search step sizes is crucial in an optimization algorithm as it determines the 
solutions diversification, exploitation performance and overall convergence process of an algorithm [20]. 
In order to better study the impact of different step size settings onto the convergence performance of the 
EM, the proposed EMs used in this experiment are also modified into two varients. These varients are set to 
operate with local search step size settings in two different extremes. EM with Larger Search Steps (EMLSS) 
is modified to search locally in a fixed search step of 0.1, while EM with Smaller Search Steps (EMSSS) is 
set to conduct its local search with a fixed search step of 0.0001. 


4. RESULTS AND ANALYSIS 

The results of the tilt angle optimization using the Ems are shown in this section. The results are 
given for the experimental PV system with the maximum output power of 210 Watts installed at the 
coordinates 3.0833° N, 101.6500° E. In order to demonstrate the impact of the solar tracking system in the 
overall power harvesting, the results are compared to that of the same system setup without any solar tracking 
system. A performance comparison with the conventional EM is included to investigate the improvements 
made by the modifications to the algorithm. 
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Table 2 shows the results comparison of the output power produced by different methods. The value 
of N is set to 100 due to the fact that the peninsular of Malaysia has the maximum solar ray around March to 
April every year. The cloudiness coefficient (F,) is set at 0.98. It can be observed from Table 2 that the 
maximum power generated by all optimization techniques is at 1200H to 1300H. The results also show that 
the output power generated by EMLSS are generally less accurate compared to all the other search 
mechanisms. Inaccurate tilt angles are kept as best solutions in EMLSS because the large search steps 
employed in the algorithm skipped some other better possible solutions. EMSSS, on the other hand, 
returned outcomes which generates higher output power compared to EMLSS. Its small search steps enabled 
it to better search the tile angle with higher accuracies. The solutions obtained by EMSSS are very 
competitive with that of the standard EM. Figure 3 gives a better illustration on the comparison of the power 
generated by each optimization technique in the form of a graph. Notice that the graph line generated by 
EMSSS is at the top of all other techniques, indicating that it generates the most power in most of the hourly 
comparisons. 


Table 2. Power Generation by Different Optimization Techniques (V=100) 








EM EMLSS EMSSS Fixed PV 
0700H 13.50W 11.03W 17.62W 3.18W 
0800H 41.84W 32.17W 39.63W 13.55W 
0900H 90.47W 87.62W 90.15W 50.78W 
1000H 133.06W 128.30W 142.99W 98.83 W 
1100H 195.63W 190.07W 192.52W 158.36W 
1200H 203.26W 202.92W 204.15W 200.98W 
1300H 199.04W 199.60W 201.90W 186.86W 
1400H 198.69W 195.34W 196.02W 163.10W 
1500H 180.36W 178.67W 183.10W 105.04W 
1600H 141.92W 135.63W 139.88W 61.77W 
1700H 101.57W 90.14W 99.73W 29.86W 
1800H 52.16W 41.96W 73.55W 12.55W 
1900H 15.30W 12.33W 12.97W 5.13W 





Results Comparison 


160 7 ~ 


0700H 0800H 0900H 1000H 1100H 1200H 1300H 1400H 1500H 1600H 1700H 1800H 1900H 


—--— EM EMLSS = -=-s+se0+ EMSSS — -- FxedPV 


Figure 3. Hourly performance comparison of the EM variants. 


Figure 4 shows the corresponding hourly tilt angles calculated by all the optimization techniques 
used in this experiment. Over-and under-compensations can be observed in the comparison of all the 
corresponding tilt angles, which in turn results in less power generated than the best optimum point. 

Figure 5 shows the convergence curves of all the search techniques in 1200H. It can be observed 
from the curves that EMLSS converges rapidly and reaches near optimal tilt angle in relatively earlier 
iterations. However, the maximum objective value returned by EMLSS is comparatively low, which indicates 
that the output power generated by EMLSS is relatively less than other techniques. EMSSS search in details 
from the beginning of the iterations. This enabled the algorithm to search deeper for a more fine-tuned tilt 
angle. However, it is clearly shown in Figure 5 that the convergence rate of EMSSS is much slower 
compared to all other benchmark techniques. 
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Tilt Angle Comparison 
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Figure 4. Comparison on the hourly tilt angles returned by the conventional EM, EMLSS, and 
EMSSS 


Some convergence processes are sampled and analyzed. The convergence performance comparison 
in Figure 5 shows that EMLSS progressed very rapidly in early stages and found near-optimal values in 
relatively earlier iteration than all other optimization techniques. However, it is found that the near optimum 
results are less accurate compared to that of the conventional EM and EMSSS. EMSSS, on the other hand, 
converge slowly but steadily towards the optimal tilt angle search. The generated power returned by the 
EMSSS is relatively higher than the EMLSS, indicating that the EMSSS has the ability to hit solutions with 
higher accuracies. 


Convergence Process 


——eEM EMLSS — - - EMSSS 


Figure 5. Convergence process of all techniques in 1200H 


5. CONCLUSION 

This paper proposes a new artificial intelligent approach to optimize the tilt angle of a solar tracking 
system using EM. In order to investigate the effect of search step size settings on the accuracy and overall 
performance of the algorithm, the EM is modified into EMLSS and EMSSS. The experimental results 
indicate that the proposed procedure performed well in searching for the best optimum tilt angles. EMSSS 
outperformed other variants in terms of accuracy and overall convergence performance, which in turn 
generates higher output energy. We thus conclude that the propose method works well in optimizing the tilt 
angle of a solar tracking system. Application of the EM in maximum power point tracking (MPPT) will be 
considered in the future work. 
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